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ABSTRACT 


Simplified expressions for the parameter controlling the load point compliance and 
strain energy release rate were obtained for the Edge Crack Torsion (ECT) specimen 
for mode III interlaminar fracture toughness. Data reduction methods for mode III 
toughness based on the present analysis are proposed. The effect of the transverse 
shear modulus, G 2 3, on mode III interlaminar fracture toughness characterization was 
evaluated. Parameters influenced by the transverse shear modulus were identified. 
Analytical results indicate that a higher value of G 2 3 results in a lower load point 
compliance and lower mode III toughness estimation. The effect of G 2 3 on the mode 
III toughness using the ECT specimen is negligible when an appropriate initial 
delamination length is chosen. A conservative estimation of the mode III toughness 
can be obtained by assuming G 2 Jl =G l2 for any initial delamination length. 


KEYWORDS: composites, compliance, strain energy release rate, mode III fracture 
toughness. 


Nomenclature 

a 

A ij 

b 

C 

d 

Aj 


delamination length, m 
extensional stiffness coefficients, N/m 
specimen width, m 
loading point compliance, m/N 
loading arm, m 

bending and twisting stiffness coefficients, Nm 



fi.fi 

Gij 

Gj 

G„, 

Gmc 


h 

Hn,Hio, In, ho 

l 


P 

Pc 

Qn 


So 


S 1 


x,y, z 


Young's moduli, GPa 
functions of (Aij, Dij, b and a ), m 
shear moduli, GPa 

total strain energy release rate, kJ/m 2 
mode III strain energy release rate, kJ/m 2 
critical mode III strain energy release rate (mode III toughness), 
kJ/m 2 

sublaminate thickness (half the specimen thickness), m 

functions of (Aij, Doe, b and a), m 

effective specimen length, m 

number of ±45 and +45 pairs in the sublaminate 

transverse load, N 

transverse load at delamination growth, N 
3-D transformed reduced stiffness, N/m 2 

characteristic root for undelaminated region, 1/m 
characteristic root for delaminated region, 1/m 
Cartesian coordinates 


Greek Letters (Symbols) 

otij parameters that relate the end loading moments to bending and 
twisting curvatures, Nm 2 
5 load point displacement, m 

8 C load point displacement associated with delamination growth, m 
0 twisting angle per unit length, or twist, 1/m 
0 fiber angle with respect to x axis for a ply in the laminate 
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INTRODUCTION 


Delamination in composites is generally characterized in terms of three fracture 
modes: mode I (opening), mode II (in-plane shear) and mode III (anti-plane shear). 
Analysis and test techniques for delamination in mode I and mode II have been 
studied extensively. O'Brien and Martin [1] summarized the mode I standardization 
efforts. Mode II fracture toughness characterization [2,3] Jiave_also received 
extensive attention, and standardization development is under way. However, 
relatively little work has been reported for mode III delamination testing. A mode III 
test is needed for complete characterization of the three fracture modes to completely 
describe mixed mode fracture processes. 

Several test methods have been proposed for measuring mode III toughness. 
Donaldson proposed a split cantilever beam method for mode III delamination test 
[4]. Becht and Gillespie used a crack rail shear test for mode III interlaminar fracture 
toughness characterization [5]. Li and Armanios developed a shear deformation 
model for the analysis of torsion loaded unidirectional and cross-ply laminates with 
double mid-plane free edge delaminations, and the total strain energy release rate was 
found to be pure mode III [6]. Shear deformation is important not only because the 
shear modulus (G^) °f composites is relatively small in comparison with in-plane 
modulus but also because the out-of-plane warping due to torsion loading is 
predominate in rectangular cross-sectioned laminates. 

Recently, Lee proposed an Edge Crack Torsion (ECT) method for mode III 
toughness characterization [7]. The test is based on a [90/(±45) n /(+45) n /90] s laminate 
specimen, with one mid-plane free edge delamination, subjected to torsion. The 
torsional response and fracture toughness of the ECT specimen were analyzed on 
the basis of the Classical Lamination Theory (CLT). However, two outstanding issues 
associated with the ECT test must be addressed. First, the presence of off-axis plies 
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will induce twist-bending coupling, and consequently, the mode II contribution to 
the total strain energy release rate needs to be identified. Second, the large difference 
(16% and 40% for n - 3 and 4, respectively) between torsional stiffness predicted by 
the CLT and that obtained from experiment [7] raises the need for a more rigorous 
analytical treatment. To address these issues, Li and Wang developed a shear 
deformation theory to analyze the ECT specimen [8]. The twist-bending coupling 
effects on the torsional stiffness and the total strain energy release rate were isolated. 
It was concluded that the mode II contribution to the total strain energy release for 
the [90/(±45) n /(+45) n /90] s class of laminates is negligible. 

In the theory developed in reference 8, in addition to the in-plane elastic properties 
of a typical unidirectional ply (laminae), the transverse shear modulus G 2 3 is needed 
to determine G UI of the ECT specimen. However, the transverse shear modulus is 
much more complicated to determine than the in-plane properties, such as E u , E 22 , 
and G 12 . Both theoretical and experimental efforts have been devoted to the 
determination of G 2 3 . Hashin and Rosen developed equations to calculate an upper 
and a lower bound of G 23 based on fiber and matrix properties [9]. Kurtz and Sun 
tested thick laminates in torsion to determine the transverse shear modulus [10]. 
However, due to lack of complete fiber and matrix properties [9] and additional 
efforts of special specimen preparation and testing [10], the transverse shear modulus 
is not easily determined for most composite materials. 

In this paper, the effect of variations in G 23 on the shear deformation theory 
prediction for the compliance and the mode III fracture toughness are investigated. 
By taking advantage of the special lay-ups for mode III testing, the lengthy 
expressions for compliance and strain energy release rate in reference 8 are simplified. 
The simplified expressions are in terms of parameters with well-known physical 
meanings. In addition, this paper also proposes data reduction methods using a single 
specimen. 
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MODE in TOUGHNESS 


A schematic of the ECT specimen is shown in Figure 1 . The delamination is located at 
the midplane on one edge and runs through the entire length of the laminate. The 
specimen is supported at three points and is loaded by a transverse load ( P ) at one 
point as shown in Figure 1 . The loading arm and the effective length of the laminate 
are denoted by d and /, respectively. This loading- condition -introduces-twisting to 
the specimen, and is approximated by the torsion loading shown in Figure 2. 
According to the St.-Venant principle, the solution of Figure 2 can be applied with 
satisfactory accuracy in region at some distance from the loading points on the 
specimen shown in Figure 1 [11]. The effect of local irregularity in stress distribution 
around loading points is neglected in the present analysis. This assumption appears 
to be justifiable because as shown in reference 8, the predicted load point compliance 
compared favorably to the experimental results. A three dimensional finite element 
analysis may be needed to determine the transition from the actual deformation and 
stress state near the loading points to the pure torsion state assumed in Figure 2. 


The laminate in Figure 2 is under a generalized plane deformation [12], and a closed 
form solution for the problem in Figure 2 has been developed in reference 8 and is 
based on a shear deformation theory. The torsional response and total strain energy 
release rate are expressed as [8] 
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where the twisting moment M T in Figure 2 is equal to Pd, (Xu a = 9oCn/3tf, etc., 0 is 
the twisting angle per unit length, or the twist, and G T denotes the total strain energy 
release rate. The parameters a i} in Equations 1 and 2 are functions of the sublaminate 
stiffness coefficients (A ;j , B^, D l} ) and specimen configurational parameters ( a , b, h) 
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and are given in reference 8. For a valid mode III toughness test, the mode II 
contribution must be eliminated, or at least minimized. The bending-twisting coupling 
should be small for a desirable lay-up (i.e., a 21 and a 21 a should be negligible). For 
such lay-ups, only the first term in the right side of Equations 1 and 2 need to be 
considered, and equations 1 and 2 may be simplified as follows 


M 


T _ 


e 

G 


= 8a 


n 


49 


1IL ~ —CL 
2 “11 ,a 


(3) 

(4) 


In addition, the sublaminates above and below the delamination plane should be 
symmetric about their own midplanes to eliminate any residual thermal stress 
contributions to strain energy release rate for delamination growth. For the symmetric 
sublaminate lay-ups with negligible bending-twisting coupling, the parameter x 
can be simplified as [8] 


“ll “(As6 + V^66 - + ^66 f\ + (^66 + 'T^66)/2 (5) 

where h is the sublaminate thickness (see Figure 2). In Equation 5, the first two terms 
represent the Classical Lamination Theory (CLT) solution, and the third and fourth 
terms contain functions f x and f 2 that represent the out-of-plane shear deformation 
contributions. These functions are given by 


h = (/„ + /,0 )(l - e- 5 ° lb - a) ) 2 + -l) 


where 


H n = D 66{*-[ 1 + J tanh ( i o( i - fl ))] 4 ' Sl ‘ I }/ A 
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For a sublaminate of thickness h, the stiffness coefficients in the above equations are 
obtained from 


(Ay **)* 

where Qy is the 3-D transformed reduced stiffness as defined in reference [13]. 


(15) 


The inverse of the loading point compliance, 1/C, is defined as the division of the load 
point displacement, 8=0W, by the transverse load, P, and can be written using 
Equations 3 and 5 as 

^ = ^ = j^2{(D 66 +^-A 66 )b-^A 66 a + D 66 f l+ (D 66+ ^A 66 )f 2 } (16) 


The mode III toughness, G IIIC can be obtained by solving Equation 3 for 0 and 
substituting into Equation 4 as 

a n.q 

(17) 


G///C = “ 


< P 2 c d ^ 
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where P c denotes the applied transverse load corresponding to delamination growth. 


TRANSVERSE SHEAR MODULUS G 23 EFFECT 

The in-plane material properties and geometry of the ECT specimen used in reference 
7 are presented in Table 1. The effect of G 2 3 on the compliance (C) and toughness 
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(G ///c ) are reflected through the parameter a u and its derivative with respect to crack 
length a. The stiffness coefficients which depend on G 23 are A u , A 55 and A 45 as seen 
from their (CLT) definitions. 


A44 - 

I-(1/2)[ G 23 COs2 + G 13 Sinl 4>]* 

(18) 

a* 

<-* 

II 

i-(i/2) [ G 23 Sin2 $ + G 13 COS 2 4 ))dz 

(19) 

11 

Tf 

1-(1/2)( G 13 - ^23 ) cos <}) sin §dz 

(20) 


where <}> denotes the angle between the fiber direction in a ply and the x axis. 

These coefficients only appear in the expressions of the characteristic roots, j 0 and s x 
as seen in Equations 13 and 14. To investigate the G 23 effects on the predictions of 
the inverse of the compliance and mode III toughness, the [90/(±45) n /(+45) n /90] s 
laminates with n = 3 and 4 were selected. Experimental results [7] showed that for the 
[90/(±45) n /(+45) n /90] s laminates, with n=l (12 plies) and 2 (20 plies), the load versus 
deflection behavior was highly nonlinear, but with n=3 (28 plies) and 4 (36 plies), 
the load versus deflection behavior was linear up to the point of delamination 
growth. These lay-ups have negligible Mode II contributions as demonstrated in 
reference 8. To assess the influence of G 23 , two extreme values of G 23 was examined; 
G 23 = 0 and G 23 =G 12 . In both cases, v 23 =v 12 is assumed to isolate the effect of G 23 
alone, even though the transverse isotropic condition, v 23 =£ 2 2 /(2G 23 )-l, is violated. 

The theoretical predictions of the inverse of the load point compliance, 1/C, for both 
the G 23 = 0 and G 23 =G 12 cases are plotted in Figures 3 and 4 for n= 3 and 4, 
respectively, along with the test data from reference 7. For both lay-ups, n = 3 and 4, 
G 23 = 0 yields a lower value for 1/C, while G 23 =G 12 gives a higher value for 1/C. The 

mode III toughness predictions and the results from the Compliance Calibration (CC) 
method used in reference 7 are given in Table 2. The delamination length and critical 
load for each lay-up used in Table 2 were taken from actual tested ECT specimens in 
reference 7. In contrast to the inverse of the load point compliance, a higher value for 
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the mode III toughness is predicted when G 23 =0 is used than when G 23 =G 12 is used. 
The lower values of the predicted mode III toughness show good agreement with the 
CC results [7] for both n = 3 and 4. Also appearing in Table 2 is a more realistic 
situation [10] where G 23 is 60% of G n and v 23 is calculated from the transverse 
isotropic condition. The calculated toughness value for this case is slightly higher 
than that corresponding to G 23 =G 12 . In conclusion, the approximation of G 23 =G 12 
appears to give a good conservative prediction of the mode III toughness. Hence, 
accurate determination of G 23 may not be mandatory to characterize the mode III 

interlaminar fracture toughness using the ECT specimen. 

COMPLIANCE CALIBRATION AND DATA REDUCTION METHODS 


The strain energy release rate can also be calculated from the compliance of the 
delaminated specimen [14] as 


r _p 2 zc 

G/// " 21 da 


( 21 ) 


If the strain energy release rate is established experimentally from measuring the 
compliance of the cracked body, the method is usually called the compliance 
calibration method. Two procedures are often adopted in the experimental 
evaluation to approximate the derivative of C with respective to a. In the first 
procedure, the delamination is extended by a small increment A a and the change in 
compliance AC is measured. Then the ratio AC/Aa is used to approximate the 
derivative in Equation 21. In the second procedure, the compliance C as a function 
of crack length a is established experimentally. A polynomial curve fit is applied to 
the experimental data. The compliance expressed in this way can then by substituted 
into Equation 21. 


For the present specimen shown in Figure 1 , it is difficult to generate a compliance 
curve from a single ECT specimen due to its unstable crack growth. A large number 


9 



of ECT specimens with different crack lengths need to be tested to obtain the 
compliance curve. Although, the shear deformation theory discussed above gives an 
effective prediction of the mode III toughness, it is not suited for data reduction 
purpose as it requires lamina elastic properties. However, it is possible to incorporate 
the shear deformation theory into the compliance calibration method for the ECT 
specimen. This may be achieved as follows. 


Substitute the derivative of a n with respect to a (using Equation 5) into Equation 4 
to obtain 



Equation 22 is plotted for the entire range of a/b (from 0 to 1 ) in Figures 5 and 6 for 
the [90/(±45) n /(+45) n /90] s class of laminates with n - 3 and 4. These results are 
plotted for three values of G 23 ; G 23 =G 12 , G 23 =0.6G 12 , G 23 = 0 in the figures. The 
dotted horizontal line in these figures represents the delamination size independent 
term (h 2 A 66 /4) of Equation 22. For a wide range of a/b (0.2 <a/b< 0.5), the parameter, 
Gjj/( 40 2 ), is nearly independent of the crack length. Furthermore for the case when 
G 23 =Gi 2 is assumed, there appears to be no variation when 0.2<a/b<0.5. Therefore 
for delamination sizes within this range, the mode III strain energy release rate may be 
approximated by 



(23) 


where, <5 C is the measured load point displacement associated with delamination 
growth, and A 66 represents the average shear modulus of the sublaminate times its 
thickness. If A 66 can be measured conveniently, the mode III toughness may be 
determined from Equation 23. 
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CONCLUDING REMARKS 


Compliance and strain energy release rate are expressed in terms of well-known 
physical parameters for desirable mode III test lay-ups for the Edge Crack Torsion 
(ECT) test. The influence of the transverse shear modulus G 23 is introduced through 
the characteristic roots, s t and s 0 , to the compliance and strain energy release rate of 
the ECT specimen. Analytical investigation ..of -the— mode -III lay-up 
[90/(±45) n /(+45) n /90] s with n= 3 and 4 indicates that higher value of G 23 results in 
lower load point compliance and lower mode III toughness estimation. In addition, 
the lower values of the predicted mode III toughness show good agreement with the 
compliance calibration results given in the literature. By assuming G 23 =G 12 , the 
present analysis gives a conservative estimation of the mode III toughness, G///, for 
any initial delamination length. Furthermore, the parameter, G ni , normalized by the 
square of the twisting angle per unit length, (G lt /Q 2 ), is nearly independent of the 
crack length when the crack length-to- width ratio, ( a/b ), is between 0.2 and 0.5. An 
alternative to the compliance calibration technique is proposed based on the 
delamination size independent term in the parameter, Gjj/(4Q 2 ). 

Acknowledgment This work was performed while the first author was a National 
Research Council Research Associate at NASA Langley Research Center. 
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Table 1. Elastic properties and configurational parameters of Carbon/Epoxy 

composite [7] 


Properties 

Configurational Parameters 

E n = 165 GPa 

b=3S.l mm (Width) 

£22 = £33 = 1 0.3 GPa 

d= 31 .8 mm (Moment arm) 

G 12 = G 13 =5.5 GPa 

1=762 mm (Effective length) 

V 12 = V 13 =0.28 

1=0.13 mm (Ply thickness) 


fr=Half laminate thickness 


Table 2. Comparison of mode III fracture toughness for 
[90/( ±4 5 ) n /( +45 ) n /90] s lay-ups. 



G,i, C (kJ/m 2 ) 

CC [7] 

Shear Deformation Theory 

n a/b 

P C (N) 


6 ^ 23=0 Gi3=G\2 G 23 =0.6G 12 

3 0.353 

943 

1.14 ±0.22 

1.42 1.18 1.23 

4 0.507 

934 

0.93 ±0.08 

1.12 0.89 0.95 
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3.0 1 


G 23 = G 12 

□ Experimental Data [7] 



Figure 3. Inverse of compliance as a function of normalized crack length a/b for 

[90/(±45) n /(+45) n /90] s lay-up with n=3. 
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Figure 4. Inverse of compliance as a function of normalized crack length a/b for 

[90/(+45) n /(+45) n /90] s lay-up with n=4. 
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